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Abstract. Planet-disk interaction predicts a change in the orbital elements of an embedded 
planet. Through linear and fully hydrodynamical studies it has been found that migration is 
typically directed inwards. Hence, this migration process gives natural explanation for the pres- 
ence of the 'hot' planets orbiting close to the parent star, and it plays a mayor role in explaining 
the formation of resonant planetary systems. 

However, standard migration models for locally isothermal disks indicate a too rapid inward 
migration for small mass planets, and a large number of massive planets are found very far 
away from the star. Recent studies, including more complete disk physics, have opened up new 
paths to slow down or even reverse migration. The new findings on migration are discussed and 
f^) ' connected to the observational properties of planetary systems. 
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1. Introduction 



CM 

O i_ The orbital elements of the observed extrasolar planets are distinctly different from 
the solar system. The major solar system planets have a nearly coplanar configuration 
and orbits with small eccentricity. In contrast, the exoplanet population displays large 
■ eccentricities, and many planets orbit their host star on very tight orbits. Recently, it has 
been discovered that inclined and retrograde orbits are quite frequent as well, at least 
for close-in planets. Historically, it was exactly the relatively 'calm' dynamical structure 

t— I | of the solar system that led to the hypothesis that planets form in protoplanetary disks. 

Within the framework of the sequential accretion scenario planet formation proceeds 
along a series of substeps, growing from small dust particles all the way to the giant 
gaseous planets. The discovery of the hot planets, which could not have been formed 

^vq i in-situ due to the hot temperatures and limited mass reservoir, gave rise to the explo- 
ration of dynamical processes that are able to change the location of planets in the disk, 
«) | accompanying the regular formation process. 

_^ ■ In the context of moons embedded in the ring system of Saturn it had been noted 

t-H | that disk-satellite interaction can alter the orbital elements of the perturbing moon, in 



particular its semi- major axis ( Goldreich fc Tremaine (1980)| . Later it was recognized 



that via the very same process, operating between a young embedded planet and the 
protoplanetary disk, it is possible to bring a planet that has formed at large distances 



from the star to its close proximity (Ward(1986) I. Hence, this migration process has 



provided a natural explanation for the population of hot planets, and their mere existence 
has been considered as evidence for the migration process. At the same time it was 
noted that planet-disk interaction may lead to eccentricity as well as inclination damping 



(Ward(1988) Ward & Hahn(1994)) 



Another indication for a planetary migration process comes from the high fraction 
(nearly 20%) of configurations in a low order mean-motion resonance, within the whole 
sample of multi-planet systems. As the direct formation of such systems seems unlikely, 
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only a dissipative process that changes the energy (semi-major axis) is able to bring 
planets from their initial non-resonant configuration into resonance. Since resonant cap- 
ture excites the planetary eccentricities typically to large values in contradiction to the 
archetypical system GJ 876, it has been inferred that planet-disk interaction should lead 
to eccentricity damping (Lee & Peale (2002) I. 

However, population synthesis models strongly indicated that the standard migration 
scenario yields very rapid inward migration rates that appear to be in disagreement with 
the observations. Additionally, recent observations of close-in planets on eccentric and 
inclined orbits have questioned the general validity of the migration paradigm to form 
them. In this review, I will first explain the basic mechanism of migration, present new 
findings on the migration rate, and then discuss its applicability with respect to the 
overall planet formation process. 



2. Origin of migration 

An embedded object disturbs the ambient disk dynamically in two important ways: 
First it divides the disk into an inner and outer disk separated by a coorbital (horse shoe) 
region. Secondly, the propagating sound waves that are sheared out by the Keplerian 
differential rotation generate density waves in the form of spiral arms in the disk. The 
created structures in the perturbed coorbital region and in the spiral arms back-react 
on the planet and cause a change in its semi-major axis. Thus, physically speaking, 
planetary migration is caused by the effect of spiral arms and corotation region. Let us 
discuss these effects in turn. 

Spiral arms: To put it simple, the spiral arms can be considered as density enhance- 
ments in the disk that 'pull' gravitationally on the planet. This gives rise to so called 
Lindblad torques that change the planet's angular momentum. For circular orbits the 
disk torque exerted on the planet is directly a measure of the speed and direction of 
migration. The inner spiral forms a leading wave that causes a positive torque, while 
the outer wave generates a negative contribution. The combined effects of both spirals 
determine then the sign and magnitude of the total torque. A positive total torque will 
add angular momentum to the planet and cause outward migration. On the other hand, 
a negative torque will induce inward migration. It turns out that under typical physical 
disk conditions the contributions of the inner and outer spiral arm are comparable in 
magnitude. However, the effect of the outer spiral quite generally wins over the inner one 
causing the planet to migrate inward. 

Corotation region: As viewed in the corotating frame, material within the corotation 
region performs so called horseshoe orbits. Here, the gas particles upon approaching the 
planet at the two ends of the horseshoe are periodically shifted from an orbit with a semi- 
major axis slightly larger than the planetary one to an orbit with slightly smaller value, 
and vice versa. Hence, at each close approach with the planet there is an exchange of an- 
gular momentum between (coorbital) disk material and the planet. The total corotation 
torque is then obtained by adding the contributions from both ends of the horseshoe. 
To obtain a net, non-zero torque requires non- vanishing radial gradients of vortensity 



and entropy across the corotation region (Baruteau & Masset(2008) I. For an ideal gas 



without friction or heat diffusion mixing effects within the horseshoe tend to flatten out 
these gradients yielding a vanishing corotation torque, or so called torque saturation. 

2.1. Type- 1 migration 

Small mass planets do not alter the global disk structure significantly, in particular they 
do not open a gap within disk. Hence, the combined effect of Lindblad and corotation 
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torques can be calculated for small planetary masses using a linear analysis. The out- 
come of such linear, no-gap studies has been termed type-I planet migration. Due to the 
complexity of considering heat generation and transport in disks these studies have relied 
nearly exclusively on simplified, locally isothermal disk models. Here, the temperature 
is assumed to be independent of height and is given by a pre-described function of ra- 
dius, T = T(r). Typically, it is assumed that the relative scale height H of the disk is a 
constant, H/r = const., yielding T cx r -1 . The total torque T to t is given as the sum of 
Lindblad and corotation torque T to t = Fl+Tcr- The speed of the induced linear, type-I 
migration scales inversely with the disk temperature (i.e. disk thickness) as oc (H/r)~ 2 , 
linear with the planet mass cx m p , and with the disk mass oc to<j. Linear models have 
been calculated for flat 2D disks as well as full 3D configurations. The problem of 2D 
simulations lies in taking into account approximately the neglected vertical stratification 
of the disk, which is typically done through a smoothing of the gravitational potential 
near the planet. Additional problems arise when considering radial gradients. Hence, 2D 
and 3D results may well yield agreeing migration rates at a particular radius but oppo- 
site dependence on radial gradients. Additionally, non-linear effects may set in already 
at a planetary mass of about 10 earth masses. New full 3D, nested grid locally isother- 
mal hydrodynamic simulations of planet-disk interaction give very good agreement with 



previous 3D linear results (Tanaka & Ward (2004)) and yield the following form for the 



total torque for small mass planets below about 10 M Earth (D'Angelo & Lubow (2010) I 



r tot = -(1.36 + 0.62a E + 0.43a T ) fe) (y) V p r*S%. (2.1) 

In eq. (|2.ip the index p refers to the planet, as and «t refer to the radial variation of 
density and temperature, such that S(r) oc r _QS and T(r) oc r~ QT . 

2.2. Type-II migration 

For larger planet masses, a gap is opened in the disk because the planet transfers angular 
momentum to the disk, positive exterior and negative interior to the planet. The depth 
of the gap that the planet carves out depends for given disk physics (temperature and 
viscosity) only on the mass of the planet (Fig. [TJ left panel). Because the density in the 
coorbital region is reduced, the corotation torques are strongly affected and are no longer 
of any importance for larger planet masses. For very large masses even the Lindblad 
torques are reduced yielding a slowing down of the planet (Fig. Q] right panel). This 
non-linear regime has been coined the type-II regime of planetary migration; here the 
drift of the planet is dominated by the disk's viscous evolution. The dip in the migration 
rate at around m p = lOMEarth has been discovered by Masset et al. (2006)| and can be 



attributed to a change in the flow structure in the vicinity of the planet due to non-linear 
effects. 

2.3. Migration in radiative disks 

The migration rates obtained through linear analysis as well as fully non-linear hydro- 
dynamical models have resulted in approximate formulae for the migration speed a of 
a planet (as quoted above for small mass regime, eq. 12. ip that are frequently used in 
population synthesis models, i.e. growth models of planets that include disk evolution 
and planet migration in parameterized form as well. These synthesis models can be used 
to calculate the final location of many planets in the mass/semi-major axis diagram 
and compare the results statistically with the observed distribution. Through variation 
of individual parameter of the model, their relative importance can be estimated. The 
results indicate in particular, that the migration for small mass planets in the type-I 
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Figure 1. Left: Azimuthally averaged surface density profile of the disk for a variety of planet 
masses, quoted in solar masses. The density structure is obtained for an isothermal disk with 
H/r = 0.05 using a constant viscosity, equivalent to a a — 0.004 at the radius of the planet, 
as obtained with 2D hydrodynamic simulations. Right: The migration rate, quoted in terms 
of the migration time scale t — a/a for different planet masses. The results refer to 3D nested 
grid sim ulations where the sym bols denote different grid layouts. The red solid line denotes the 
result of |Tanaka et al. (2002) 



regime is by far too fast to account for the observed distribution, the majority of plan- 



ets would have been lost to the star, as shown for example by Ida & Lin (2008) and 



Mordasini et al. (2009) Only a significant reduction in the type-I migration speed gives 
satisfactory results. Suggested remedies included: stochastic migration of a planet in a 
turbulent disk, migration in inviscid self-gravitating disks or nonlinear effects. 

Here, we shall concentrate on a very simple and straight forward improvement of the 
models, that represent a possible solution to the type-I migration problem: the inclusion 
of more accurate physics. As mentioned above, past modeling relied nearly exclusively 
on the simplified locally isothermal models, which has the advantage that no energy 
equation has to be considered. Taking more realistic thermodynamics into account re- 
quires the incorporation of a heating and cooling mechanism. The importance of radiative 
diffusion has first be pointed out by Paardekooper & Mellema (2006) and recent pa- 



pers quote approximate formulas for viscous and diffusive disks (Masset & Casoli(2010) 



Paardekooper et al.(2010) I. To demonstrate the effect for two-dimensional flat disks, we 



show results for a planet-disk simulation where we include viscous heating, local radiative 
cooling as well diffusive radiative transport in the disk's plane. The energy equation then 
reads 

dY, ^ T + V ■ (£c v Tu) = -pV u +D-Q - 2HV ■ F (2.2) 
ot 

where £ is the surface density, T the midplane temperature, p the pressure, D the 
viscous dissipation, Q the radiative cooling and F the radiative flux in the midplane. 
Models where the various contributions on the right hand side of eq. (|2.2p were selectively 



switched off and on, have been constructed by Kley & Crida (2008) 



The effect of this procedure on the resulting torque is shown in the left panel of Figj2j 
The basis for all the models is the same equilibrium disk model constructed using all terms 
on the rhs. of eq. (I2.2[) and no planet. Embedding a planet of 20M eart h yields in the long 
run positive torques only for the radiative disks, where the maximum effect is given when 
only viscous heating and local radiative cooling are considered. The inclusion of diffusion 
in the disk midplane yield a slightly reduced torque. Since the initial state consists of 
a non- vanishing negative radial entropy gradient, the adiabatic model shows a positive 
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Figure 2. Time evolution of the torque acting on a 20 M ear th planet embedded in a disk with 
0.01 solar masses with a radial range of 0.4 to 2.5 rj up . Left: Simulations using a constant value 
of the viscosity. From bottom to top the curves indicate simulations where i) no energy equation 
(isothermal), n) only the first term on the rhs. of eq. (|2.2[) (adiabatic), Hi) all terms on the rhs. 
(fully radiative), and iv) all but the last term on the rhs. (heating/cooling), have been used, 
(after Kley & Crida, 2008) Right: The influence of viscosity on the resulting torque for fully 
radiative simulations using an alpha type viscosity. In the simulations the dissipation has been 
kept fixed, varying only the value of a in the momentum equation. 

torque during the first 20 orbits directly after insertion of the planet. However, in the long 
run the torque becomes negative as in the isothermal case because the material within the 
horseshoe region is mixed thoroughly, wiping out the entropy gradient. Hence, the positive 
corotation torque is obliviated and the negative Lindblad contribution dominates. The 
adiabatic case does not approach the isothermal result because the corresponding sounds 
speeds are different. The right hand panel of Fig |5] demonstrates that a non- vanishing 
value of the viscosity is necessary to maintain torque desaturation. These models have 
used an a-type viscosity in contrast to the results displayed on the left side. To compare 
results directly, it has been ensured that the thermal state of the models is identical 
despite the different value of a used in the angular momentum equation. 
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Figure 3. Left: Radial torque density, T(r), for fully 3D radiation-hydrodynamical simulations 
with an embedded 20M ear th planet. The data are shown after equilibrium has been reached. 
Different assumptions for the thermodynamical state are made, (after Kley et al., 2009) Right: 
The total torque acting on the planet for a 2D fully radiative simulation for a variety of planet 
masses, (after Kley & Crida, 2008) 



2.4. Fully three-dimensional radiative disks 

The previous 2D cases have been repeated using full 3D radiative disk simulations with 
an identical physical setup. These clearly confirmed the existence of outward migration in 
the presence of radiative diffusion. To compare isothermal, adiabatic and fully radiative 
simulations, all start from the same initial state which corresponds to an equilibrium 
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of the fully radiative case including viscous heating and radiative diffusion. The spatial 
origin of the torques can be analyzed using for example the radial torque density T(r) 
which is defined through T tot = jT(r)dr, where T tot is the total torque acting on the 
planet. A plot of the radial torque distribution (Fig. [31 left panel) shows that at t = 80 
the corotation torques have saturated in the isothermal and adiabatic case, and only 
the Lindblad contributions remain. Obviously, the net effect is the sum of two contri- 
butions that have opposite sign and are of comparable magnitude. The negative part of 
the outer spiral arm has a slightly larger amplitude than the inner contribution. Again, 
the isothermal and adiabatic torques differ due to the different sound speed. The fully 
radiative case agrees for radii larger than a p with the adiabatic model while there ex- 
ists a well pronounced torque maximum just inside of the planet. This contribution is 
responsible for the torque reversal. The right panel of Fig. [3] shows that the strength 
of this positive corotation effect also scales with the square of the planet mass up to 
about 20 to 25 M eart h- Beyond this mass, gap opening begins and only the Lindblad 
torques remain, and above 40 M eart h planets begin to migrate inwards again. In the full 
3D simulations the results are qualitatively the same, the averaging procedure, that is 
necessary in 2D, leads to some quantitative differences Kley et al. (2009) Interestingly, 
the full 3D results show even a stronger effect. New population synthesis models based 
on the modified migration rates indicate better agreement with the observational data 
set (eg. Mordasini, this volume). 



3. Eccentricity and inclination 

In addition to a change in semi-major axis, planet-disk interaction will modify the plan- 
etary eccentricity (e) and inclination (i) as well. Extending previous isothermal studies by 
Cresswell et al. (2007)[ fully 3D radiative disk simulations have been performed recently 



to study the evolution e and i. The results indicate that both are damped for all plane- 
tary masses ( Bitsch fc Kley (2010)| and this volume). For small values of e and i that are 
below about 2H/r the damping occurs exponentially on timescales comparable to the 
linear estimates by Tanaka & Ward (2004) For larger values, damping is slowed down 
and follows approximately e oc e~ 2 , and for the damping of i an identical relation holds, 
surprisingly. Interestingly, the presence of outward migration is coupled to the magni- 
tude of e and i. Outward migration only occurs for eccentricities smaller than about 0.02, 
and inclinations below about 4°. This reduction is due to the fact that for non-circular 
orbits the flow structure in the corotation region becomes strongly time dependent and 
no stationary corotation torque can develop. More details about the evolution of inclined 
and eccentric planets in 3D disks are given in |Bitsch fc Kley (2010)[ and Bitsch & Kley 
(this volume). 



4. Resonant systems 

The mere existence of resonant planetary systems is a strong indication that dissi- 
pative mechanisms changing the semi-major axis of planets must have operated, as the 
likelihood of forming planets in these configurations in situ or later through scattering 
processes seems to be small. On the other hand, convergent differential migration of a 
pair of planets will lead under quite general conditions to capture into resonance. The 
presence of an inner disk in combination with realistic inner boundaries leads to very 
good agreement with the observations for the best observed system GJ 876, as demon- 
strated in Crida et al. (2008) and shown in Fig.UJ Through transit timing measurements 
the first double eclipsing system, Kepler 9, has been discovered, and the data clearly 
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Figure 4. Evolution of the semi-major axis a and eccentricity e of a pair of planets with physical 
parameter resembling the observed planetary system GJ 876. For the initial 75 years the planets 
are kept fixed and thereafter they are allowed to migrate freely within the disk. The index 1 
refers to the inner and 2 to the outer planet. The horizontal lines refer to the observed values of 
GJ 876. The evoluti on has been stopped after the planets have reached the observed distances 
from the star, (after Crida et al. (2008) I 



indicate a low order mean motion resonance, probably 2:1 (Holman et al.(2010)|. The 
proximity of the planets to the star and the near coplanarity of the system strongly hints 
towards a migration scenario for the formation of the system. 

The importance of dynamical migration models for systems of planets is indicated by 
the system HD 45364, where two planets engaged in a 3:2 resonance have been discovered 



by Correia et al.(2009) The inferred orbital parameters for the two planets are semi- 
major axes of ai = 0.681AU and ai = 0.897AU, and eccentricities of e\ — 0.168 and 
ei = 0.097, respectively. Fully non-linear hydrodynamical planet-disk models have been 



constructed for this system by Rein et al. (2010) For suitable disk parameter, the planets 



enter indeed into the 3:2 resonance through a convergent migration process. After the 
planets have reached their observed semi-major axis, a theoretical RV-curve has been 
calculated. Surprisingly, even though the hydrodynamically obtained eccentricities (ei = 
0.036, e 2 = 0.017) are quite different, the model fits the observed data point equally 
well as the published data, see Fig. [5] and Rein et al. (2010) This pronounced dynamical 
difference between the two models that match the observations equally well, can only 
be resolved with more observational data. Hence, the system HD 45364 is a very good 
example that it is urgently necessary to have sufficient and better observational data for 
interacting multi-planet systems to obtain the orbital parameters, that are required to 
constrain the theoretical models. 



5. Summary 

We have shown that just by inclusion of more accurate physics, namely radiation trans- 
port, it is possible to reduce significantly the otherwise too rapid type I migration. The 
effect is driven by corotation torques, and the requirement to maintain the horseshoe 
torques unsaturated is the presence of viscosity and radiative cooling (or diffusion) , with 
timescales of the order of the libration time of the horseshoe material near the separatrix. 
Including this effect, population synthesis models yield better agreement with the obser- 
vations. We find that eccentricity and inclination are typically damped by planet-disk 
interaction for all planet masses. 

Differential migration in multi-planet systems frequently result in resonant capture. 
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Figure 5. Synthetic radial velocity curves HD 45364 together with observational data points. 
The light green curve is the fit as given in the discovery paper (|Rein et al. (2010)] ) . The darker 
blue curve has been obtained by performing full hydrodynamical simulation of the system. Both 
curves have comparable \ 2 values. 

Upon capture the eccentricity of the planets is strongly increased. For continued migra- 
tion the systems remain stable only when eccentricity is damped by the disk. We have 
shown that the formation of the systems GJ 876 and HD 45364 can naturally be ex- 
plained by planet-disk interaction migration scenarios. For resonant systems where the 
disk action is reduced in the final stages of the planet formation process the left over 
configuration may result in an unstable system. The last phase will then be dominated 
by scattering processes which may pump up the planetary eccentricities and possibly 
their inclinations to the large values observed. 
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